Tunable transmission and bistability in left-handed bandgap structures 
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We study the defect-induced nonlinear transmission of a periodic structure created by alternating 
slabs of two materials with positive and negative refractive index. We demonstrate bistable switching 
and tunable nonlinear transmission in a novel type of bandgap that corresponds to the vanishing 
average refractive index, and compare the observed effects for two types of the bandgaps. 



Materials with both negative electric permittivity and 
negative magnetic permeability were suggested theoreti- 
cally a longtime ago and they were termed left-handed 
materials |2|. Such materials can also be described by a 
negative refractive index, as was demonstrated by several 
reliable experiments 0, |j] and numerical finite-difference 
time-domain simulations 0. 

Multilayered structures composed of left-handed (LH) 
materials can be considered as a sequence of flat lenses 
that provide an optical cancellation of the conventional 
right-handed (RH) layers leading to either enhanced or 
suppressed transmission 0, Q • More importantly, a one- 
dimensional stack of layers with alternating RH and LH 
material with vanishing average refractive index <n> 
displays a novel bandgap 0, El thai differs from a 
conventional Bragg reflection gap. 

In this letter, we study defect-induced bistable switch- 
ing in the novel bandgap structures and, for the first time 
to our knowledge, compare the effects observed for zero- 
index and Bragg reflection gaps. The particular struc- 
ture we study by numerical pseudospectral time-domain 
(PSTD) simulations consists of seven periods of a LH-RH 
double-layer (indicated in Fig.^). Each individual layer 
has equal width a. The LH layer of the central period is 
replaced by a Kerr-type nonlinear material, which con- 
stitutes a structural defect of the periodic system even in 
the linear regime. We also study the effects produced by 
a change of the defect position for two types of bandgaps. 

We calculate the electric and magnetic fields directly 
from Maxwell's equations using the PSTD method llll. 
The differential Maxwell equations are approximated by 
difference equations and solved iteratively. The material 
properties are treated through the electric permittivity 
s r and magnetic permeability /i r . The PSTD method is 
advantageous for the modeling of interfaces where both e r 
and fj, r change because it samples all material properties 
at the same location [l^ . 

The linear properties of the LH material are described 
by Lorentz dispersion characteristics in both e r and /i r , 
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FIG. 1: Layout of the system with snapshots of the elec- 
tric field amplitude. A Gaussian pulse with carrier frequen- 
cies lo c = 2tt x 14.93 GHz (solid) and lo c = 2?r x 25.88 GHz 
(dashed) is incident on the structure. The fields are shown at 
the pulse peak. The shaded background indicates the mate- 
rial structure, with dark grey denoting the Kerr material and 
light grey the LH slabs. The defect layer is centered within 
the structure. The dotted vertical line shows the source loca- 
tion, the dash-dotted vertical lines denote the boundaries of 
the considered structure. The field is incident from the left. 



D(w) = e r (cj)£oE(w), B(w) = /j. r (u))fioH(w), trans- 
formed into ordinary differential equations in the time 
domain, and subsequently approximated by difference 
equations |l| that are incorporated into the PSTD al- 
gorithm. This allows one to model frequency dependent 
material properties in PSTD simulations. The nonlinear 
Kerr material is described by a dielectric function 
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These functions are substituted into the relations, 



There are several established methods of including the 
nonlinear material response in the PSTD algorithm [hH 
Ell 113 Here, we use the electric field of the previous 
time step to evaluate Eq. © 15, 16]. 

We also calculate the properties of the structure using 
the transfer- matrix method (TMM) 18]. This method 
allows an exact analytical solution of the linear problem. 
In particular, one can relate incident, transmitted and 
reflected fields using the transfer matrix of the structure, 
and thus obtain an explicit expression for the transmis- 
sion and reflection coefficients. 

We describe the LH material by Eqs. CQ,© with pa- 
rameters chosen to give a refractive index n ~ — 1 at 
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FIG. 2: Linear amplitude transmission coefficient of the struc- 
ture calculated by the TMM (line) and PSTD (symbols) 
methods. The band gaps are indicated by shaded regions. 
The <n> — (m = 0) bandgap lies around 15 GHz. The 
first conventional Bragg gap (m = +l)lies around 25 GHz 
while the first Bragg gap in the left-handed regime (m = — 1) 
lies around 9 GHz. 

the design frequency f = 15 GHz. We use ui pe = 
1.1543 x 10 11 s~\ Loie = io Xm = 2ir x 5 x 10 6 s -1 , 
and ui pm = 1.6324 x 10 11 s _1 , and include small losses, 
7 e = 2 x 7 ? „ = 2n x 6 x 10 6 s -1 . With these parame- 
ters the LH meta-material is left-handed for frequencies 
/ < 18.5 GHz and right-handed for / > 26 GHz. The 
slab thickness a is chosen to be Ao/4, where Ao is the free- 
space wavelength at f . We use air as the RH medium. 
The parameters for the Kerr material are X = 3 and 
•^(3) = 4 m 2 /v 2 . In each case the PSTD simulations use 
a discretization of 100 points per Ao and a time step cor- 
responding to half the Courant stability limit 11] of the 
linear case. 

To validate our PSTD calculations, we first study the 
amplitude transmission spectrum in the linear regime, 
shown in Fig. [3 and find excellent agreement between 
the TMM and PSTD methods. For frequencies below 
8 GHz, the discretization level in the PSTD simulations 
is insufficient for accurate results, due to the high refrac- 
tive index of the LH material, and some discrepancies 
appear. The structure exhibits band gaps due to Bragg 
scattering, both in the RH and in the LH frequency re- 
gion. An additional band gap appears around the fre- 
quency where the average refractive index vanishes, as 
was addressed earlier 0, 0, E| • These band gaps can be 
identified through the index to in the usual Bragg condi- 
tion, 

fcRHORH + fcLHOLH = WWT, (4) 

where kj is the respective wave number and dj is the 
respective layer thickness. The band gap index to can 
be any integer, including zero and negative numbers. We 
find that the defect layer introduces a transmission peak 
into each of the shown band gaps. The remainder of this 
letter will focus on the comparison of the <n> — gap 
(to = 0) and the first RH Bragg gap (m = +1). 
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FIG. 3: Dependence of the defect frequency on the position 
offset A of the defect layer in the linear limit. 
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FIG. 4: Time-integrated transmitted intensity per pulse, rel- 
ative to that of the incident pulse, in the m = gap (solid) 
and the m = +1 gap (dashed). Bistability occurs in the m — 
gap for all offsets, while the bistability region for the m = +1 
gap is indicated in the graph by grey shading. 

We calculate the frequency of the defect modes in the 
band gaps for different locations of the defect layer in 
our structure in the linear shown in Fig. [3J The 

defect layer is shifted by a quantity a A (|A| < 1) from 
the central position (A — 0) between its two neighboring 
LH layers without changing the rest of the structure. It 
is remarkable that in the to — gap the defect frequency 
does not depend on the defect position while the one in 
the to = +1 gap is shifted with changing A. When the 
thickness of the defect layer is not equal to a (not shown) , 
the to = defect frequency also depends on A. This 
dependence, however, is much weaker than that of the 
to = +1 defect. 

In the nonlinear regime, the behavior of the fields de- 
pends on the intensity of the electric field inside the 
nonlinear layer. We study this problem by PSTD sim- 
ulations. The incident field has a Gaussian envelope in 
time with a width parameter of 1528//o and amplitude 
0.2 V/m. The carrier frequency is lo c — 2ir x 14.93 GHz 
for the m — band gap and ui c = 2tt x 25.88 GHz for the 
to = +1 gap. These frequencies lie below their respective 
defect frequencies at A = with equal relative detuning. 

In Fig. 01 we show the relative time-integrated trans- 
mitted intensity per pulse versus the position of the non- 
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FIG. 5: Output vs. input intensity of the structure, showing 
hysteresis behavior for defect offsets of A = 0, —0.4 in the 
m = gap and A = 0, -0.08, -0.4 in the m = +1 gap. The 
A = —0.4 curve in the m = +1 gap is an almost vanishing 
straight line without bistability. 



linear layer in the linear and nonlinear regimes for pulses 
at both carrier frequencies. The two bands display very 
different character. In the linear regime, the transmission 
in the m = gap is only weakly dependent on A with 
a maximum at A = 0. In the m = +1 gap, the linear 
transmission shows a double peak structure near A = 0. 
Away from this region the transmission is much reduced. 
For A rs ±1 the transmission increases again. In the 
nonlinear case the transmission is generally reduced. In 
the m = gap, it maintains its character with a single 
peak in the center, which actually exceeds the peak of 
the linear case. In the m = +1 gap, the double peak 
character is lost and a much smaller single peak in the 
center and weak wings for large layer offsets remain. For 
all other A the transmission almost vanishes. 

One of the most promising applications of periodic 
structures with embedded nonlinear defect layers is the 
possibility to achieve intensity-dependent tunable trans- 



mission in the spectral bandgaps usually associated with 
bistability. Optical bistability is a powerful concept that 
could be explored to implement all-optical logic devices 
and switches. In nonlinear systems that display bistabil- 
ity, the output intensity is a strong nonlinear function of 
the input intensity [T^. [2(j ■ 

In the m = gap, we observe bistability for all offsets 
while in the m = +1 gap, bistability is observed for very 
small values of |A|. These effects are related to the be- 
havior of the defect frequency in the m = +1 gap with 
different A. In this gap the defect frequency shifts away 
from the carrier frequency with increased A and eventu- 
ally approaches it again. In contrast to this, the defect 
frequency in the m = gap stays very near to the carrier 
frequency for all layer offsets. In Fig. [5] we show the out- 
put vs. input intensity during the time of the pulse. The 
left panel shows the behavior in the m — gap for offsets 
A = and A = —0.4, while the right panel shows the be- 
havior in the m = +1 gap for offsets A = 0, A = —0.08, 
and A = —0.4. The bistable behavior, where the output 
intensity with increasing input intensity is different from 
the one with decreasing input intensity, is clearly visible. 
The threshold switching intensity in the m = +1 gap is 
significantly lower than in the m — gap. On the other 
hand, the difference between switch-up and switch-down 
threshold is much larger in the m = gap. The differ- 
ence between the behavior in the two gaps is illustrated 
by the A = —0.4 curves. While the m = gap shows sig- 
nificant bistability, the m = +1 gap shows no bistability 
and almost no transmission. 

In conclusion, we have studied the defect-induced 
nonlinear transmission and bistability for two types of 
bandgaps of a layered structure composed of two mate- 
rials with positive and negative refractive index with a 
Kerr nonlinear defect. We demonstrated a number of 
unique features of the bandgap associated with zero av- 
erage refractive index of the structure. 
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